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ABSTRACT 

DNA methylation has been proven to be a critical 
epigenetic mark important for various cellular 
processes. Here, we report that redox-active 
quinones, a ubiquitous class of chemicals found in 
natural products, cancer therapeutics and environ- 
ment, stimulate the conversion of 5mC to 5hmC 
in vivo, and increase 5hmC in 5751 genes in cells. 
5hmC increase is associated with significantly 
altered gene expression of 3414 genes. 
Interestingly, in quinone-treated cells, labile iron- 
sensitive protein ferritin light chain showed a signifi- 
cant increase at both mRNA and protein levels 
indicating a role of iron regulation in stimulating 
Tet-mediated 5mC oxidation. Consistently, the 
deprivation of cellular labile iron using specific 
chelator blocked the 5hmC increase, and a 
delivery of labile iron increased the 5hmC level. 
Moreover, both Tet1/Tet2 knockout and dimethylox- 
alylglycine-induced Tet inhibition diminished the 
5hmC increase. These results suggest an iron- 
regulated Tet-dependent DNA demethylation mech- 
anism mediated by redox-active biomolecules. 



INTRODUCTION 

DNA methylation regulates gene expression and renders 
cellular identity, and also is critically involved in genome 
imprinting, inactivation of X chromosome and parasite 
elements, and documentation of epigenetic memory (1). 
The 5-methylcytosine (5mC) could be oxidized by the 
ten-eleven translocation (TET) family of hydroxylases to 
the sixth DNA base 5-hydroxymethylcytosine (5hmC) 
(2-6). Varying degrees of the 5hmC modification has 
been found in several tissue types (3,7,8). Genome-wide 
mapping analysis has shown a strong enrichment of 
5hmC within exons and near transcription start sites 
(TSSs), pointing to a potential role in transcriptional regu- 
lation (9-11). Furthermore, impaired hydroxylation of 
5mC has been observed in myelodysplastic syndromes 
and several forms of myeloid leukemias (12,13). 

Currently, the generation of 5hmC in mammalian cells 
is primarily attributed to enzymatic oxidation of 5mC 
by either of the three existing TET proteins (4-6,14). 
The generated 5hmC could be further converted to remov- 
able 5-formylcytosine and 5-carboxycytosine by further 
oxidation, providing a mechanism for active DNA 
demethylation (5,6). Interestingly, small biomolecules 
can alter the epigenetic DNA methylation patterns and 
regulate life processes (14), even cause some diseases like 
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cancer (15-18). These xenobiotics include metals (19-21), 
vitamin C (14), endocrine-disrupting chemicals (22) and 
some persistent organic pollutants (23,24). On the other 
hand, the alteration of epigenetic states including DNA 
methylome could be exploited to treat cancers (25). 
However, little is known about the effects of small bio- 
molecules on Tet-mediated DNA demethylation. 

Since the oxidation of 5mC catalyzed by Tet proteins 
is assisted with some small cofactors, e.g., iron (II) and 
a-ketoglutarate (4-6), the catalytic activity of Tet 
proteins should be affected by the levels of iron (II) and 
a-ketoglutarate in mammalian cells. The cellular levels of 
iron (II) and a-ketoglutarate are altered by cellular milieu, 
in particular redox status (26-30). Quinones are a ubiqui- 
tous class of redox-active compounds found in natural 
products, cancer therapeutics, endogenous biochemical, 
and environment pollutants or generated through metab- 
olism of hydroquinones and/or catechols (31). We specu- 
late that redox-active quinones may change the hemostasis 
of iron or a-ketoglutarate, therefore influence the catalytic 
activity of Tet dioxygenases. Since Tet dioxygenses mediate 
the oxidation of 5mC and initiate passive and active DNA 
demthylation (5,6), it expects that quinones may alter 
DNA demethylation process through changing Tet oxida- 
tion activity. Here, we tested if quinones could impact the 
conversion of 5mC to 5hmC at cellular levels. If so, we 
would characterize the possible roles of those compounds 
in modulating DNA demethylation and investigate the 
underlying mechanism. 

MATERIALS AND METHODS 

Cell culture and treatment 

Human fetal lung fibroblast cell line MRC-5, human lung 
adenocarcinoma cell line A549 and hepatocellular carcin- 
oma cell line HepG2 (supplied by cell culture center, 
Institute of Basic Medical Sciences, Chinese Academy of 
Medical Sciences, Beijing, China) were cultured in 
Dulbecco's modified Eagle's medium (DMEM) with high 
glucose medium and in RPMI 1640 medium, respectively, 
in 5% C0 2 at 37°C, both of which contained 10% fetal 
bovine serum, 100 U/ml penicillin, 100 ug/ml streptomycin. 
The wild type (WT) and Tetl and Tet2 genes knocked out 
(Tetl/Tet2- / -) mouse embryo stem cell were maintained in 
DMEM medium (HyClone, Thermo Fisher Scientific Inc., 
MA, USA) supplemented with 20% ES FBS (Gibco, Life 
Technologies Corporation., Grand Island, NY, USA), 
0.1 mM non-essential amino acids, 2mM L-glutamine, 
0.1 mM P-mercaptoethanol, 1000 U/ml leukemia inhibi- 
tory factor (LIF) (Millipore, Billerica, MA, USA), 
1.0 uM PD 0325901 (Stemgent, Cambridge, MA, USA) 
and 3.0 uM CHIR 99021 (Stemgent, Cambridge, MA, 
USA). For all experiments, mES cells were trypsinized 
and plated in culture dishes pretreated with 0.1% gelatin, 
then incubated in a humidified 37°C incubator supplied 
with 5% C0 2 . For the generation of 5hmC in cells, 
MRC-5 or mES cells were treated with 20 uM or 50 uM 
tetrachloro-l,4-benzoquinone (TCBQ) or tetrachloro-1,4- 
hydroquinone (TCHQ) for 24 h. For the effect of 
dimethyloxaloylglycine (DMOG) on TCBQ-dependent 



generation of 5hmC, A549 cells were treated by 50 uM 
TCBQ or 200 uM DMOG or both for 24 h. For the dose- 
dependent assay, A549 cells (5 x 10 5 cells) were seeded in 
culture medium for 24 h and then treated with 0, 1, 5, 20 
and 50 uM of TCBQ or TCHQ for 24 h. The medium was 
removed and the cells were washed three times with phos- 
phate buffered saline (PBS), and fresh medium containing 
the same concentration of TCBQ or TCHQ added. The 
treatment was consecutively repeated three times over a 
total period of 72 h. To determine the effect of exposure 
time on the level of 5hmC, A549 cells were treated with 
50 uM TCBQ or TCHQ for 0, 6, 12, 24 and 48 h (two 
24 h treatments), 72 h (three 24 h treatments). In order to 
elucidate the repair mechanism of 5hmC in A549 cells, after 
being treated with 50 uM TCBQ for 72 h as described 
above, the TCBQ-containing medium was removed and 
the cells were washed three times with PBS. TCBQ- 
treated A549 cells were then recultured in fresh medium 
without chemical for 0, 12, 24, 48, 72 and 96 h. HepG2 
cells were treated with TCBQ (20 uM) and TCHQ 
(50 uM) for 72 h and the treatment was carried out as 
described above. In all, 50 uM of tetrabromo-1,4- 
benzoquinone (TBrBQ), 2,5-dichloro-l,4-benzoquinone 
(2,5-DCBQ), 2,6-dichloro-l,4-benzoquinone (2,6-DCBQ), 
2-chloro-l,4-benzoquinone (2-CBQ), O-chloranil and 2- 
methyl-l,4-benzoquinone (2-MBQ) were also used in 
A549 cells exposure experiments for 24 h. All the cultured 
cells treated with chemicals were harvested for further 
analysis, and the exposure experiments were repeated 
three times. 

DNA extraction and enzymatic digestion 

Genomic DNA was extracted from the harvested cells 
using a Genomic DNA Purification Kit (Promega, 
Madison, WI, USA) according to the manufacturer's 
instructions. To prevent artificial oxidative damage to 
DNA, 0.1 mM DFO was added in the processes of 
DNA extraction and purification (32). The concentration 
and quality of the extracted DNA were evaluated by 
measuring the absorbance at 260 nm and 280 nm. 

The DNA (5 ug) was digested with 1 U DNase I, 2 U 
calf intestinal phosphatase and 0.005 U snake venom 
phosphodiesterase I (New England Biolabs, Ipswich, 
MA, USA) at 37°C for 24 h. A total of 1 mM DFO was 
used during the enzymatic digestion of DNA (32). 
Proteins from the DNA digestion system were removed 
prior to the UHPLC-MS/MS analysis, using Microcon 
centrifugal filter device (Millipore, Bedford, MA, USA) 
with the 3000 D cutoff membrane by centrifuging at 
10 000 g for lh. 

UHPLC-MRM MS/MS analysis 

The Agilent 1200 Series Rapid Resolution LC system and a 
reverse-phase Zorbax SB-C18 2.1 x 100 mm column 
(1.8 (im particles) were applied in the UHPLC analysis. 
The method was same as recently described (14) or 
with a minor modification. In brief, the digested DNA 
(5.0-15.0 ul) was injected onto the column and nucleosides 
separation was completed using the mobile phase of 95% 
water (containing 0. 1 % formic acid) and 5.0% methanol at 
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a flow rate of 0.3 ml/min, or a gradient elution was used for 
separation: 0-3 min, 5.0% B; 3-6 min, 15.0% B; 6-10 min, 
100% B; 10-15 min, 5.0% B. Solvent A was an aqueous 
solution of 2.0 mM NH4HCO3 (pH 9.0), and solvent B was 
100% methanol, the flow-rate was 0.25 ml/min. Mass spec- 
trometric detection was achieved by the Agilent (Santa 
Clara, CA, USA) 641 0B Triple Quadrupole mass spec- 
trometer with an electrospray ionization source. Positive 
ionization mode, 3500 V capillary voltage as well as 
300°C and 91/min nitrogen drying gas were selected for 
the experiment. For nucleotides analysis, the fragmentor 
voltage and collision energy were performed at 90 V and 
5eV, respectively, with a scan time of 100 ms. Multiple 
reaction monitoring (MRM) mode was used for the LC- 
MS/MS analysis. 5hmC in cellular DNA was detected in 
the form of the mononucleoside 5hmdC by monitoring the 
transitions of mjz 258.1 — ► 142.1. The amount of 5hmC 
was calibrated by standard curve. 

Dot blot for 5hmC analysis 

Briefly, DNA samples were diluted with lOmM Tris-HCl, 
pH 8.0, and denatured by heating at 95°C for 10 min, and 
chilled on ice for 5 min. Then the DNA samples were 
loaded on Amersham Hybond N+ membrane (GE 
Healthcare) using a 96-well dot-blot apparatus (Bio- 
Rad). After being baked in 80°C for 2h and blocked by 
5% non-fat milk for 1 h at room temperature, the 
membrane was incubated with a polyclonal anti-5-hmC 
antibody (Active Motif 39791, 1:10,000) at 4°C overnight. 
5-hmC was visualized by using chemiluminescence after 
being incubated with secondary antibody against rabbit 
(Roch, 1: 8000) for 30 min at room temperature. 5hmC 
standard (active motif) at 1, 0.5, 0.25 ng from top to 
bottom was used as positive control. 

Western blot analysis 

Cytoplasmic and nucleic extracts from MRC-5 cells 
treated by 20 uM TCBQ or 0.02% (v/v) DMSO were 
prepared with NE-PER® Nuclear and Cytoplasmic 
Extraction Reagents (Thermo Scientific) according to 
the manufacturer's instructions. Protein concentrations 
were measured using the Bradford Assay, and 60 (ig 
protein extracts were separated on a 12% sodium 
dodecyl sulfate-polyacrylamide gel electrophoresis. Then 
proteins were transferred to a nitrocellulose membrane, 
blocked with 5% non-fat milk and incubated overnight 
at 4°C with antibodies against the light chain of ferritin 
(FTL) (abeam, 1:500), ferroportin (Fpn) (sigma, 1:1000), 
GAPDH (Sigma, 1:1000). Blots were incubated with 
secondary antibody against mouse (1:8000) and rabbit 
(1: 8000) for lh before visualization. Cytoplasmic and 
nucleic extracts from MRC-5 cells treated with 100 uM 
ferric nitrilotriacetate (Fe-NTA) were applied for 
positive control. 

5hmC pull-down, deep sequencing and data analysis 

Genomic DNA from TCBQ or Control DMSO treated 
MRC-5 cells was purified by using Wizard genomic DNA 
purification kit (Promega) according to the manufac- 
turer's instructions. Genomic DNA samples were 



further sonicated into short fragments (about 300 bp) 
by using Covaris DNA shearing (Covaris S2 SonoLAB 
Single) with microTUBEs according to the manufac- 
turer's instructions. Sonicated DNA was then purified 
with phenol/chloroform/isoamyl alcohol (25:24:1) pre- 
cipitation. 5hmC was pulled down as described previ- 
ously (11). After purification, 10 ng of the pull-downed 
5hmC containing DNA was end-repaired, adenylated, 
ligated to adapters and single-end sequenced on a 
HiSeq 2000 system (Illumina) according to the manufac- 
turer's recommendations for Illumina ChlP-Seq to 
identify peak enrichment. To the 101-bp sequence 
reads, as a subset of reads may contain all or part of 
the adaptor oligonucleotide sequence, a custom pair- 
wise alignment script was applied to find reads hit by 
adaptor sequence, which required at least five bases 
match from the beginning of the adaptor, and at least 
90% similarity. If detected, the read was trimmed to the 
preceding base. Low-quality base trimming was per- 
formed by Trimmomatic (http://www.usadellab. 
org/cms/index.php?page = trimmomatic) with default par- 
ameters. Clean reads were mapped to the Homo sapiens 
reference genome (Hml9) using BWA (33) with default 
parameters. The maximum number of mismatches 
allowed was 3. To the reads mapped to different 
genome locations with the same quality of mapping, we 
randomly chose one location. The 5hmC peaks were 
analyzed using MACS (34). By default, MACS retains 
no more than one read on the same starting position in 
the genome, which can removes redundant reads. Peak 
enrichment with a linear normalization were implemented 
in MACS. Genes with fold change over 5 were selected 
for GO and KEGG analysis. The metagene analysis of 
the 5hmC peaks were performed by CEAS package 
(33,34) (http://liulab.dfci.harvard.edu/CEAS/). 

RNA-seq 

Total RNA was isolated from TCBQ treated and DMSO 
treated MRC-5 cells using RNAzol (Molecular Research 
Center, Inc.). Poly(A) RNA from 1 ug of total RNA was 
used to generate the cDNA library according to TruSeq 
RNA Sample Prep Kit protocol which was then be 
sequenced using the Illumina system. Since we sequenced 
single-end 101 bp reads, we kept a caution avoiding 
removing duplicates for that may introduce biases in the 
expression analysis. We aligned those high-quality reads 
to the masked Homo sapiens reference genome (hgl9) and 
junction sequences using BWA 0.5.8 software (35) 
allowing up to four mismatches. All the uniquely 
mapped reads were used to define gene expression level 
based on our developed method (36). We modified the 
RPKM (measured in reads per kilobase of exon per 
million mapped sequence reads, which is a normalized 
measure of exonic read density) calculation described by 
Mortazavi et al. (37) to include junction fragments. By 
including reads that map to junction fragments, we 
ensure that genes with many introns are not underrepre- 
sented in the RPKM estimate. To examine differential 
expression, we used a MARS (MA-plot-based method 
with Random Sampling model) method to identify 
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differentially expressed genes which was described in 
DNAseq (38). To determine the DEGs, the P-value and 
FDR-adjusted P-value (using both Benjamini-Hochberg 
correction (39) and the method of Storey and Tibshirani 
(40)) were both considered: P < 0.001 and q < 0.05 (indeed 
both q-values < 0.003). 

Gene ontology and enrichment analysis 

Gene ontology (GO) analysis of differentially expressed 
genes was performed using DAVID (http://david.abcc. 
ncifcrf.gov/). An enrichment map of 5hmC enriched 
genes or overlapped genes from 5hmC-seq and differen- 
tially expressed genes was generated by using the 
Cytoscape (Ver. 2.8.2) installed with Enrichment 
Map plugin. Red node represents each enriched GO 
pathway. Node size equals the total number of genes 
in each GO pathway. Edge thickness is in proportion to 
the number of overlapped genes between nodes. GO 
pathways of related functions are sorted into one cluster, 
marked with circles and labels. Gene numbers in each 
cluster are also labeled. 

Verification of 5hmC-enriched regions with boronic 
acid method 

Genomic DNA purified from TCBQ or DMSO treated 
MRC-5 cells was glucosylated with the EpiMark™ 
5hmC and 5mC Analysis Kit (NEB, Ipswich, MA, 
USA) according to the manufacturer's instructions with 
minor modification. For 5hmC-enriched regions verifica- 
tion, genomic DNA with glucosylation or not was treated 
by 5mM boronic acid (BA) reagent in lOOmM Na 2 HP0 4 
(pH 8.5) and 50 mM NaCl buffer, then the treated DNA 
was diluted to lOng/ul and 1 ul was used in the 25 ul RT- 
qPCR reaction containing lxGoTaq® SYBR Green 
qPCR Master Mix (Promega, Madison, WI, USA) and 
0.2 uM forward and reverse primers. qPCR was per- 
formed on a Mx3500P Real-time PCR system 
(Stratagene) and the cycling condition was as follows: 
95°C for 2min and then 40 cycles of PCR at 95°C for 
15s, 57°C for 15s and 25°C for 60s. The experiments 
were independently performed triplicate. The primer 
pairs used for detection of 5hmC-enriched regions were 
listed in Supplementary Table S4. 

RNA isolation, cDNA synthesis and quantitative 
real-time PCR 

MRC-5 cells were treated with 20 uM TCBQ dissolved in 
DMSO solution or DMSO solution only. Cells were 
trypsinized routinely after 24 h treatment. Total RNA 
was extracted using TRI® Reagent (Sigma) (Invitrogen). 
cDNA was synthesized by RevertAid™ First Strand 
cDNA Synthesis Kit with Oligo dT primers (Fermentas, 
K1622) following the manufacturer's recommendations. 
PCRs were carried out on a 7500 Real-time PCR 
system (Applied Biosystems) in 25 ul reaction volume 
containing 1 ul cDNA, 2X SYBR Green PCR Master 
Mix (TAKARA) and 200 nM primers. All samples were 
analyzed by RT-qPCR for three times. The primer pairs 
used for detection of transcripts were listed in 
Supplementary Table S5. 



Statistical analysis 

The f-test using two-way ANOVA in Grouped Analyses 
of Prism5 software was applied for statistical analysis. 

RESULTS 

The conversion of 5mC to 5hmC by redox-active quinones 
in vivo 

To test whether redox-active quinones could influence 
the formation of 5hmC in vivo, we treated MRC-5 cells 
with TCBQ of 20 uM. TCBQ, as one example of 
halogenated quinoid compounds, is a reactive metabolite 
of the widely used biocide pentachlorophenol (PCP) 
(41 ,42), which has been classified as a Group 2B compounds 
(possible carcinogen to human) by IARC (43). Moreover, 
TCBQ analogs have been recently identified as toxic disin- 
fection byproducts in drinking water (44). By coupling 
ultra-high performance liquid chromatography-triple 
quadrupole mass spectrometry with multiple-reaction 
monitoring (UHPLC-MRM-QQQ), we were able to 
detect trace amounts of cellular 5hmC (Supplementary 
Figure SI). Indeed, TCBQ treatment increased the level of 
5hmC by ~1. 9-fold compared to the DMSO control (from 
25.8 per million C to 48.8 per million C) (Figure 1A and B). 
Furthermore, TCHQ (20 uM), the reduced form of TCBQ, 
could also induce 5hmC formation (Figure 1A and B). An 
anti-5hmC antibody-based dot-blot analysis further 
confirmed the increase of 5hmC induced by 20 uM 
TCBQ- or TCHQ-treatment while the uploaded amounts 
of genomic DNA were same as indicated by methyl blue 
staining (Figure 1C). The increase of 5hmC observed after 
TCBQ- or TCHQ-treatment was both dose (0-50 uM) and 
time (0-72 h) dependent (Supplementary Figure S2A-D). 
The level of 5hmC in TCBQ- or TCHQ-treated cells 
increases over 3 days of exposure, indicating that the 
generated 5hmC is stable and accumulates in cells. The 
generated 5hmC by 72-h TCBQ exposure could not be 
removed within 96 h after the TCBQ was completely 
withdrawn from the culture (Supplementary Figure S2E), 
supporting the accumulation and persistence of TCBQ- 
induced 5hmC. The TCBQ/TCHQ-induced formation of 
5hmC was also detected in human A549 cells and HepG2 
cells (Supplementary Figure S2A, S2C and S2F). 

TCBQ and TCHQ tested here are reactive metabolites of 
PCP, which has been found in body fluids (plasma or urine) 
of non-occupationally exposed individuals (0-7.0 uM) 
and in the blood of occupationally exposed workers 
(22.5-170 uM) (42,45). Since >20% of PCP could be con- 
verted to TCHQ and TCBQ (46), the total concentration of 
TCBQ and TCHQ is estimated about 0-1.5 uM for non- 
occupationally exposed individuals and 4.5-34 uM for 
occupationally exposed individuals. Here we demonstrated 
that the cells treated by TCHQ at 1.0 uM could signifi- 
cantly increase the level of 5hmC (Supplementary 
Figure S2C). The tested doses of TCBQ and TCHQ here 
are comparable to that of non-occupational exposure 
and are even lower than that for occupational exposure. 
These results clearly suggest that the observed effects of 
TCBQ and TCHQ on cellular processes by 5hmC are 
physiologically relevant. 
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Figure 1. The conversion of 5mC to 5hmC by redox-active quinones at cellular level. (A and B) UHPLC-MRM MS/MS analysis (A) and quan- 
titation (B) of 5hmC in the enzymatic digest of genomic DNA from MRC-5 cells. The cells were treated with 20 uM TCBQ or TCHQ. The peak 
corresponding to 5hmC is denoted by the arrow, and the other two peaks are un-identified. (C) Immunodotblotting analysis of 5hmC in MRC-5 cells 
upon 20 uM TCBQ- or TCHQ-treatment (left panel). The amount of uploaded DNA samples was indicated by methyl blue staining (right panel). 
N.C and P.C. represent negative and positive control, respectively. Control genomic DNA was extracted from the cells only treated with the solvent 
DMSO (0.02%, v/v) with the same amount as that used in TCBQ/TCHQ treatment. **/ > <0.01, determined using Student's /-test, two-tailed. Error 
bars, mean ± S.E.M . for triplicate experiments. 



The formation of 5hmC stimulated by quinone is Tet 
dependent 

Since quinones are redox active (31), TCBQ may stimulate 
the formation of highly reactive oxygen species (ROS). 
Indeed, a significant increase in highly reactive ROS was 
observed in TCBQ-treated cells as measured by 
2',7'-dihydroxychlorofluorescein diacetate (DCFH-DA) 
fluorescence method (47) (Supplementary Figure S3). 
The production of highly reactive ROS was dependent 
on TCBQ concentration (Supplementary Figure S3C), 
and the treatment with 50 uM TCBQ increased the 
average fluorescent intensity 25-fold above control level 
after 2h (Supplementary Figure S3C). It is not known 
whether the elevated ROS can directly react with 
genomic DNA as a type of chemical reaction (48,49) or 
indirectly stimulate the enzymatically catalytic reaction 
enhancing 5hmC formation. So, we first tested whether 
the generation of 5hmC is pertinent to TET family 
dioxygenases, which are enzymatically associated with 
5mC oxidation and resultant 5hmC formation (2-6). We 



measured 5hmC in the Tet-proficient mouse embryonic 
stem (ES) cells treated with TCBQ or TCHQ. Without 
TCBQ/TCHQ treatment, the frequency of 5hmC in 
genomic DNA of ES cells is estimated as 7.39 x 10 2 per 
million C, which is about 30 times higher than that of 
MRC-5 cells. This indicates a much higher oxidation 
activity of Tet proteins in mouse ES cells. Treatment 
with 20 uM TCBQ or TCHQ over 24 h increased the 
level of 5hmC to 1.57 x 10 3 and 1.75 x 10 3 per million C, 
respectively (Figure 2A and B). Evidently, the generation 
of 5hmC stimulated by quinones increases with an 
increased Tet oxidation activity, suggesting a linkage of 
quinone-stimulation effect with Tet oxidation activity. By 
double knockout of Tetl/Tet2 (Tetl/Tet2~ ! ~) in mouse ES 
cells, 5hmC was barely detectable even after TCBQ- or 
TCHQ-treatment (Figure 2A and B). These results 
clearly support that the formation of 5hmC by quinones 
requires Tetl or Tet2 proteins, rather than a direct 
reaction of highly reactive ROS with genomic DNA. 
Moreover, dimethyloxalylglycine (DMOG), a cell- 
permeable competitor of a-ketoglutarate (50), inhibited 
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Figure 2. TET-dependent generation of 5hmC by redox-active quinones in cells. (A and B) UHPLC-MRM MS/MS analysis (A) and quantitation (B) 
of 5hmC in the enzymatic digest of genomic DNA from wild type (WT) and TetljTet2 double knock-out (TetljTet2~'~) mouse ES cells treated by 
20 uM TCBQ or TCHQ. (C and D) UHPLC-MRM MS/MS analysis (C) and quantitation (D) of 5hmC in genomic DNA from 50 uM TCBQ-treated 
A549 cells in the absence or presence of 200 uM DMOG. Control DNA was extracted from cells treated with 0.02% (v/v) DMSO. *P<0.05, 
**.P<0.01, determined using Student's (-test, two-tailed. Error bars, mean ± S.E.M. for triplicate experiments. 



the 5hmC formation in TCBQ-treated MRC-5 cells 
(Figure 2C and D). Together, the above results consist- 
ently support the involvement of Tet proteins in quinones- 
stimulated generation of 5hmC in vivo. 

The quinone stimulation of Tet oxidation activity is 
associated with iron regulation in vivo 

The next question we asked is how quinones stimulate 
the 5mC oxidation activity of Tet proteins. Since certain 
ROS could elevate the labile iron as reported previously 
(26-29), we speculated that ROS induced by TCBQ treat- 
ment might stimulate the catalytic activity of Tet proteins 



by increasing the availability of its cofactor iron (II) at 
cellular level. Interestingly, we observed an elevation of 
the light chain of ferritin (a major iron-storage protein 
at the cellular and organismal level) in TCBQ-treated 
MRC-5 at both levels of mRNA and protein 
(Figures 3A and 7C). Essentially, both the levels of cyto- 
plasmic and nuclear ferritin light chain (FTL) were 
up-regulated upon TCBQ treatment. The level of FTL is 
strongly regulated by the level of labile iron (51,52), and 
an increase in FTL is an indicator of an elevation of labile 
iron (53). Meanwhile, a cell surface iron transporter 
ferroportin (Fpn) facilitating the export of cytosolic iron 
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Figure 3. The quinone stimulation of Tet oxidation activity is associated with iron regulation in vivo. (A) Analysis of expression of ferritin light-chain 
(FTL) and ferroportin (Fpn) proteins within cytoplasmic and nucleic extracts of 20 uM TCBQ-treated MRC-5 cells. The extracts of 100 uM ferric 
nitrilotriacetate (Fe-NTA)-treated cells were applied for positive control (P.C). Control protein was extracted from 0.02% (v/v) DMSO-treated cells. 
(B) Quantifications of 5hmC in genomic DNA from MRC-5 cells treated by 100 uM ammonium ferric citrate (FAC) for 12 h. Control genomic DNA 
was extracted from the cells treated with 0.1% (v/v) sterilized ddH 2 0. (C and D) Quantifications of 5hmC in genomic DNA from A549 cells treated 
with 50 uM TCBQ (C) or TCHQ (D) in the absence or presence of 1,10-phenanthroline (0-500 uM). */ > <0.05, **P<0.01, ***P<0.001, ns / > >0.05, 
determined using Student's /-test, two-tailed. Error bars, mean ± S.E.M. for triplicate experiments. 



(53) was almost not influenced upon TCBQ treatment. 
These results indicated a higher level of cellular labile 
iron contributing to the stimulation of Tet oxidation 
activity by quinones. The increase in labile iron pool 
would supply more iron (II) to iron (II)- and a- 
ketoglutarate-dependent Tet proteins, by which the Tet- 
catalyzed reaction (5mC oxidation and 5hmC formation) 
could be enhanced. This further prompted us to investi- 
gate the role of iron availability in Tet-catalyzed reactions 
at cellular level. 

We first delivered iron (II) to MRC-5 cells with an 
assistance of a weak chelator ammonium ferric citrate 
As a result of 12 h iron delivery, the level of 5hmC effi- 
ciently increased by 30% (Figure 3B). We then examined 
the role played by intracellular labile iron in the formation 
of 5hmC using 1,10-phenanthroline, a lipophilic and 



membrane-permeable iron-specific chelating agent (54). 
We found that 5hmC was reduced by this iron-chelating 
agent (5-500 uM) to the background level in TCBQ- or 
TCHQ-treated cells (Figure 3C and D), confirming the 
involvement of intracellular labile iron. 

Generation of 5hmC in cells by other redox-active 
quinones 

Interestingly, the generation of 5hmC in cellular DNA 
could also be promoted by a number of other quinoid 
compounds, including less halogenated quinones such as 
2-chloro-l, 4-benzoquinone (2-CBQ), 2,5-dichloro-l, 
4-benzoquinone (2,5-DCBQ) and 2,6-dichloro-l,4-benzo- 
quinone (2,6-DCBQ), and highly halogenated quinones, 
e.g., tetrabromo-1, 4-benzoquinone (TBrBQ) and 
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Figure 4. The levels of 5hmC in A549 cells treated with other redox active quinones. The tested quinones include tetrabromo-l,4-benzoquinone 
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mean±S.E.M. for triplicate experiments. 



O-chloranil (Figure 4). In addition, the methylated 
quinones (e.g., 2-methyl-l,4-benzoquinone (2-MBQ)) 
also promoted 5hmC formation. As shown in 
Figure 4A, halogenated quinones had higher activities 
for the generation of 5hmC in cells compared to non- 
halogenated quinones, TBrBQ-treatment increased the 
level of 5hmC by ~2. 8-fold compared to the control, 
and ~2. 0-fold increase was observed upon O-chloranil 
treatment. However, even for the non-halogenated 
quinone, 2-MBQ, the level of 5hmC was also raised by 
~ 1.8-fold compared to the DMSO control, which was 
comparable to the level of 5hmC in cells treated with 
TCBQ (-1.9-fold, Figure IB). This indicates that a 
number of redox-active quinoid compounds have a 
strong ability to stimulate the conversion of 5mC to 
5hmC in mammalian cells (Figure 4). 

Genome-wide profile of 5hmC induced by TCBQ 
treatment 

Quinones can be derived from many natural as well as 
synthetic chemical compounds and are therefore an in- 
evitable part of the cellular environment (31). 5mC is an 
important epigenetic mark and our finding that quinones 
are capable of promoting the conversion of 5mC to 
5hmC, thereby altering this mark, might have a wide 
range of biological consequences for the cell. To test 
the extension of this quinone-induced epigenetic 
change, we analyzed the genome-wide 5hmC distribution 
after quinone treatment by deep-sequencing of chem- 
ically labeled 5hmC (11). We found that the TCBQ- 
induced 5hmC is associated with gene-rich regions in 



MRC-5 cells (Figure 5A and B and Supplementary 
Table SI). We identified a total of 17 767 5hmC peaks 
in the genome of MRC-5 cells treated with TCBQ 
(P < 10~ , fold change >5) among which 76.9% of the 
peaks were located within either intergenic regions 
(7324) or introns (6330). The remaining 5hmC peaks 
were found in exons (1117), 5'UTR (1081) and 3'UTR 
(651), respectively (Figure 5C and Supplementary Table 
S2). The 5hmC peaks that fall into either the gene 
bodies or up-/down-stream of the genes might regulate 
gene expression and were associated with 5751 genes. In 
contrast, we only identified 6976 5hmC peaks in the 
genome of control DMSO solvent treated MRC-5 cells 
(Figure 5C and Supplementary Table S2). Strikingly, 
following TCBQ treatment we observed a significant 
enrichment of 5hmC in gene bodies as well as in prox- 
imal upstream and downstream regions relative to TSSs 
and transcription end sites (TESs) (Figure 5D). Thus, 
these analyses demonstrate that TCBQ treatment 
induces genome-wide but not loci-specific 5hmC forma- 
tion. In order to investigate the potential functions of 
5hmC, we did KEGG pathway and GO term analysis of 
the 5751 genes, which revealed that these genes belong 
to various functional groups including protein catabolic 
process, phosphorylation, cell cycle, pathways in cancer, 
cellular transport, protein complex assembly, etc. 
(Supplementary Figure S4). 

TCBQ regulates gene expression via 5hmC enrichment 

To determine if TCBQ-induced 5hmC formation results 
in changes in gene expression, we took advantage of 
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Figure 5. Genome-wide 5hmC formation by TCBQ treatment. (A) 5hmC densities in the region of chrl9:49,319, 120-49,696,264 of TCBQ treated 
(blue) and Control DMSO treated (pink) MRC-5 cell samples by 5hmC pull-down and deep sequencing. (B) 5hmC densities in the region of 
chrl:8,904,958-8,950,900 of TCBQ treated (blue) and Control DMSO treated (pink) MRC-5 cell samples by 5hmC pull-down, the red arrows indicate 
the TCBQ induced peaks. (C) 5hmC peak numbers of TCBQ treated (blue) and Control DMSO treated (pink) samples in different genomic regions 
(P < 0.005, fold change > 5). (D) Normalized 5hmC tag density distribution across the gene body. Each gene body was normalized to 0%— 100%. 
Normalized Tag density is plotted from 4kb upstream of TSSs to 4kb downstream of TESs. 



RNA-sequencing technology (RNA-Seq) to compare the 
transcriptomes of Control DMSO and TCBQ treated 
MRC-5 cells (Supplementary Table S3). We detected the 
expression of a total of 17 569 genes, with 10484 genes 
showing a significant difference in expression by using 
DEGseq software in which genes with P< 0.001 and 
q < 0.05 were defined as differentially expressed genes 
(Figure 6A and Supplementary Table S6). We next 
asked whether this difference in gene expression was 



associated with 5hmC modification of the genes. Indeed, 
3414 of the above 5751 genes enriching 5hmC show a 
significant difference in expression (Figure 6B). GO and 
enrichment analysis of these genes showed that TCBQ 
treatment resulted in 5hmC-mediated gene expression 
changes of genes belonging to several different functional 
groups, including protein catabolic process, apoptosis, 
protein localization and transport, protein complex 
assembly, RNA processing, phosphorylation, cell 
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Figure 6. TCBQ regulates gene expression via increasing 5hmC enrichment. (A) Transcriptome analysis identified a total of 17 569 gene expressed 
in TCBQ and Control DMSO treated MRC-5 cells, and 10484 of them are differentially expressed in TCBQ- compared to DMSO-treated cells 
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gene expression is significantly increased in gene bodies in TCBQ- compared to DMSO-treatment (P < 0.001, q<0.05). (C) GO and enrichment 
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TCBQ treatment compared to DMSO treatment. Red nodes represent enriched GO pathways. Clusters of functionally related gene-sets were 
generated manually. The cutoff parameters for enrichment analysis with Cytoscape software are as follows: P < 0.005, FDR q < 0.05, overlap 
cutoff > 0.5. (D) Normalized 5hmC density distribution in genes with high, medium or low expression levels. Genes were sorted by expression 
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organization, regulation of cell size and other regulatory 
pathways (Figure 6C). This suggests that quinone-depend- 
ent 5hmC generation could influence a broad range of 
cellular functions. 

Metagene 5hmC read density profiles for RefSeq tran- 
scripts revealed that gene expression level correlated with 
both intragenic and proximal enrichment of 5hmC 
(Figure 6D), consistent with a role for 5hmC in maintain- 
ing and/or promoting gene expression (1 1). Among these 
genes, about 55% of the genes enriching 5hmC were 
down-regulated. The other 1523 of the above 3414 genes 
enriching 5hmC were up-regulated and selected for 
further enrichment analysis based on DAVID GO tool 
and visualized as an enrichment map using Cytoscape 
software (Figure 7A and Supplementary Figure S5A). 
A large proportion of these genes are related to protein 
catabolic process, apoptosis signaling pathway, cell local- 
ization and transport process, RNA processing 



procedures (Figure 7A). We further validated several 
genes which have increased levels in both 5hmC and 
gene expression by conventional RT-qPCR and BA 
combined with qPCR assays, respectively (Figure 7B 
and C; Supplementary Figure S5B and 5C; 
Supplementary Tables S4 and S5). 



DISCUSSION 

Here we demonstrate that redox-active chemicals stimu- 
late oxidative conversion of 5mC to 5hmC in a TET 
dioxygenase-dependent manner. Redox-active quinones 
induce genome-wide 5hmC modifications affecting expres- 
sion of thousands genes involved in a broad range of 
cellular processes. The conversion of 5mC to 5hmC has 
been linked to epigenetic reprogramming, dynamic DNA 
demethylation and regulation of tissue-specific gene 
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Figure 7. 5hmC enrichment induced by TCBQ treatment positively correlated to gene expression. (A) GO and enrichment analysis for the 
overlapped up-regulated genes is based on the DAVID GO result. The cutoff parameters for enrichment analysis with Cytoscape software are as 
follows: P < 0.005, FDR q < 0.05, overlap cutoff > 0.5. (B) Representative examples of individual genes and their values of 5hmC and gene expression 
from TCBQ treated (blue) and DMSO treated (pink) samples, FTL gene (left panel) and Ubiquitin C (UBC) gene (right panel) are shown here. 
(C) Verification of the 5hmC-enriched regions and expression of FTL gene (left panel) and UBC gene (right panel) by BA method combined with 
qPCR and RT-qPCR, respectively. The primer-targeted regions are noted by red lines in (B). */ > <0.05, **/"<0.01, determined using Student's /-test, 
two-tailed. Error bars, mean ± S.E.M. for triplicate experiments. 



expression (55-57). Our data suggest that redox-active 
chemicals by inducing 5hmC formation could change 
the mammalian epigenome and hence affect various 
cellular processes. On the other hand, redox-active 
quinones can promote genome-wide DNA demethylation 
through enhancing genome-wide 5hmC formation as a 
result of Tet-mediated 5mC oxidation. 

Aberrant 5hmC generation by quinones might very 
likely change the normal DNA methylation landscape, 
and thus alter the epigenetic control of gene expression. 
Genome-wide hypomethylation has been observed in 
cancer cells (58) and DNA demethylation mediated 
oncogene activation has been proposed to contribute to 
tumorigenesis (59). Furthermore, hypomethylation has 
been associated with genomic instability (60). 

Interestingly, polychlorinated phenols, found in 
widely used biocides, generate redox-active halogenated 
quinones when metabolized and are under suspicion of 
being carcinogenic (41,61). We speculate that one cause 
of these biocides potential carcinogenicity could be their 
ability to convert 5mC to 5hmC in DNA. Finally, epigen- 
etic manipulation by redox-active chemicals might aid 
in the development of new cancer therapies and regenera- 
tive medicine in the future. 
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